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We describe a rigorous construction, using matched asymptotic expansions, which establishes 
under very general conditions that local terrestrial and solar-system experiments will measure the 
effects of varying ‘constants’ of Nature occurring on cosmological scales to computable precision. In 
particular, ‘constants’ driven by scalar fields will still be found to evolve in time when observed within 
virialised structures like clusters, galaxies, and planetary systems. This provides a justification for 
combining cosmological and terrestrial constraints on the possible time variation of many assumed 
‘constants’ of Nature, including the fine structure constant and Newton’s gravitation constant. 
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Investigations into the possibility that the fine struc¬ 
ture ‘constant’, cn, is slowly varying [l| combine observa¬ 
tional evidence from quasar spectra with indirect solar 
system data and laboratory constraints on allowed varia¬ 
tions. In so doing, most commentators invariably assume 
that the local and cosmological observations are directly 
comparable ;2]. However, this is a strong assumption 
that requires proof. Why should local experiments ‘see’ 
variations in ‘constants’ occurring on cosmological scales 
in regions which are gravitationally unbound and partic¬ 
ipating in the universal Hubble flow? In this letter we 
describe the first rigorous demonstration of when such 
an assumption is justified and identify the key inequality 
that must be satisfied if local variations of ‘constants’ are 
to track global variations. 

Theories which introduce varying constants self- 
consistently into Einstein’s conception of a gravitation 
theory do so by deriving the ‘constant’ C from some 
scalar field, <f>, so C —> <C(</>). The variations of this 
scalar field contribute to the spacetime curvature like 
other forms of mass-energy. They must also conserve 
energy and momentum and so their dynamics are con¬ 
strained by a wave equation 

= X ^fA4>) L k{p,p), (1) 

i,j,k 

where <j> is associated with the variation of a ‘constant’ C 
via a relation C = /(</>); the constants A, are dimension¬ 
less measures of the strength of the space-time variation 
of C, fj (cp) are functions determined by the definition 
of </>, and L k (p, p) are some linear combinations of the 
density, p, and pressure, p, of the matter that is coupled 
to the field <j>. This form includes all the standard the¬ 
ories for varying constants, like those for the Newtonian 
gravitation ’constant’ G, a, and the electron-proton mass 
ratio, described in refs. 0 |H Q ■ 

We shall refer to our scalar field as the ‘dilaton’, denote 
it by <fi, and analyse the form of CD further by assuming 
that 4>(x, t) satisfies a conservation equation that can be 
reduced to the form 

( 2 ) 


where T is the trace of the energy momentum tensor, 
T = Tjf, with any contribution from a cosmological con¬ 
stant, A, excluded; we absorb any dilaton to cosmological 
constant coupling into the definition of V(4>). The dila¬ 
ton to matter coupling, B {(/>), and the self-interaction 
potential, V ((/>), are arbitrary functions of <f>, k = 8irG 
and c = h = 1. This covers a wide range of theories 
which describe the spacetime variation of ‘constants’ of 
Nature, and includes Einstein-frame Brans-Dicke (BD), 
|3i], and all other, single field, scalar-tensor theories of 
gravity. In cosmologies that are composed of dust, cos¬ 
mological constant and radiation it will also contain the 
Bekenstein-Sandvik-Barrow-Magueijo (BSBM) of vary¬ 
ing a , [f|, and other (single dilaton) theories which de¬ 
scribe the variation of standard model couplings, Jgj. 
We assume that the background metric is of Friedmann- 
Robertson-Walker (FRW) form 

ds 2 = dt 2 — — a 2 (t)r 2 { dd 2 + sin 2 9d (j) 2 }, (3) 

1 — kr z 

where k is the curvature parameter. In the background 
universe the dilaton is also assumed to be homogeneous 
(so <f> = 4> c (t), where c denotes background quantities) 
and satisfies: 

(a 3 <p c ^) = B t< j,((f> c )K (e c - 3p c ) - (<j) c ) . (4) 

The dilaton conservation equation also reduces to an 
ordinary differential equation if we are only interested 
in its time-independent mode cf>(r) in a static spherically 
symmetric Schwarzschild spacetime. In the actual uni¬ 
verse, however, spacetimes that look static in some lo¬ 
cality must match smoothly on to the cosmological back¬ 
ground on large scales. Hence, to model the evolution 
of the dilaton accurately in some inhomogeneous region, 
embedded in a homogeneous background, we require that 
at large distances 4> —> </> c appropriately. Even if spher¬ 
ical symmetry is assumed, the conservation equation is 
difficult to solve, either exactly or approximately. Even 
numerical models are technically difficult to set-up, see 
ref. {§], and only allow us to consider one particular 
choice of B(^>), V((j)), and the spacetime background at 


□</> = B^{<j))KT - V t4> {(j)) 
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a time. If we are to bring all the experimental evidence 
to bear on these models, we need to know when we are 
allowed to use local observations to constrain cosmologi¬ 
cal variations of ‘constants’ like a. Hence, we would like 
to find the conditions which ensure that the local time- 
variation of 4> tracks the cosmological one 

<fi(x,t) & <j> c (t) (5) 

to some specified precision. 

We shall assume that the dilaton field is only weakly 
coupled to gravity, and has a negligible effect on the 
background spacetime geometry. If we ignore this back- 
reaction then we can consider the dilaton evolution on 
a fixed background with given matter density. We will 
take the background to be an exact solution to Einstein’s 
equations with matter possessing the following proper¬ 
ties. The local metric is approximately Schwarzschild, 
with mass m at time t = to, inside a closed region of 
spacetime outside a surface at r = R s . The metric for 
r < R s is left unspecified. Asymptotically, the metric 
must approach FRW and the whole spacetime should 
become exact FRW in the limit m —> 0. As the lo¬ 
cal inhomogeneous energy density of the asymptotically 
FRW spacetime tends to zero, the local spacetime met¬ 
ric exterior to r = R s must tend to the Schwarzschild 
metric. Similar scenarios have been considered by Wet- 
terich in E3 and in a very recent work by Avelino et al 
eh Whilst our conclusions will be similar to the ones 
claimed in both of these studies, we believe that neither 
study proves what is claimed. In we identified an flaw 
in 10]. When this is corrected the proof provided breaks 
down. A key point of the matched asymptotic expansion 
method is that the excluded terms in the asymptotic ex¬ 
pansions for <f> and <j> must be shown to be smaller than 
the included ones. These are non-trivial requirements. 
In E3 they are simply assumed and so no proof of the 
main claim is provided. 

In the local, approximately Schwarzschild region, the 
intrinsic interior length scale, A/, of a sphere centred on 
the Schwarzschild mass with surface area InR 2 is defined 
by the Riemann invariant: 

p 3 /2 

o = = t-Att (6) 

(2m) ' 

In the asymptotically FRW region, the intrinsic length 
scale is proportional to the inverse root of the local en¬ 
ergy density: l/^/ne c + A, where e c is the total energy 
density of matter, We shall assume that the FRW region 
is approximately flat (k = 0), and we define a length scale 
appropriate for this exterior region at epoch at t = to 
equal to the inverse Hubble parameter at that time: 

L e = 1/H 0 . 

For realistic models Le A/; we therefore define S to 
be a small parameter given by 6 = Li/Le- 


We solve the dilaton equation 0 using an asymp¬ 
totic expansion for small S in the interior region, where 
LJ 1 (t — to), Rj 1 r ~ 0(1), and also in the exterior where 
Hot ~ 0(1) ~ Hor. In general, an asymptotic expan¬ 
sion valid in one region will not also be valid in the 
other. We can only enforce a subset of our boundary 
conditions in each region, hence the resulting expansions 
will each contain undefined constants of integration. We 
have used the method of matched asymptotic expan¬ 
sions to remove this ambiguity in the expansions by as¬ 
suming that there exists an intermediate region, where 
Lj\t - t 0 ),Rj 1 r ~ 0(d~ a ) and H 0 t, H 0 r ~ 0(S 1 ~ a ), 
in which both the interior and exterior expansions for 
<f> remain valid. Then, by the uniqueness of asymptotic 
expansions, the two expressions must be equal there, see 
|12j . In this way we can determine the previous unknown 
constants of integration in each expansion, and effectively 
apply all of our boundary conditions to each expansion. 

Formally, we require that, as 5 —> 0, our choice of 
spacetime background should be FRW at zeroth order 
in the exterior expansion and Schwarzschild to lowest or¬ 
der in the interior expansion. All spherically-symmetric 
solutions to Einstein’s equations with pressureless matter 
and a A term fall into the Tolman-Bondi class (see ref. 
0). They are parametrised by two arbitrary functions 
of one spatial variable, r. We have only looked at the flat, 
Tolman-Bondi models with non-simultaneous initial sin¬ 
gularities El El , and the non-flat Tolman-Bondi models 
with simultaneous initial singularities. These two sub¬ 
cases are fully specified by prescribing only one function: 
the matter density e(r,L) on some initial space-like hy¬ 
persurface, tj. Another metric used to study the effect of 
the universe’s expansion on solar system dynamics is the 
McVittie metric [13l fTs ]. In the exterior limit the McVit- 
tie metric asymptotes to a dust-plus-A FRW cosmology. 
For radii r where Hq 1 >r> 2m, the McVittie metric 
looks like Schwarzschild; the horizon in the McVittie met¬ 
ric possesses a naked curvature singularity and so cannot 
model a black hole in an expanding universe. 

To zeroth order in both the exterior and interior, all 
these background choices will look alike and the solution 
for (j> will also be the same in all cases. As boundary con¬ 
ditions we take —> (f> c {t) as r —> oo, where </> c (t) satisfies 
the cosmological dilaton equation 0. We assume that 
the dilaton flux out of the ‘star’, with surface at R = R s , 
is given, to leading order, by 

( 2m \ — 

~ R S [ l - Sufio\ R=Rs = 2?t iF (4 > 0 ) 

dR’ R ,2 B^((j)o(R' ,t))ne(R'), 

where R is the ‘physical radius’, i.e. a surface (t, R) = 
const has surface area 47ri? 2 . For black-holes we must de¬ 
mand F (ip 0 ) = 0, we otherwise expect F (<fo) ~ B^(<j) c ). 
At higher orders we find the flux by perturbing the above 
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expression as explained in llffl . The zeroth-order exterior 
solution is then just <t> = <j> c {i)- In the interior we find 

0 (o) = 4>e {ST) - F (4>o) In ^1 - |^) , 


where T = L~J l {t — t 0 ) and £ = R~ l R. The matching 
procedure gives <j> e {ST) = (j> c {t). 

We have performed our analysis under the assumption 
that the background universe is described by McVittie 
metric. In isotropic coordinates this is 


ds 2 = 


1 - fj,{t,r) 


1 + n{t,r)_ 


+r 2 {dd 2 +sm 2 ed(j> 2 }} , 


where 


^• r) = 2^( 1 + L ‘ kr ^ 


1/2 


( 8 ) 


and m and k are arbitrary constants. In the limit a = 1, 
m is the Schwarzscliild mass, k is the curvature of the 
surfaces (t. r) = const in the m = 0 (i.e. FRW) limit. 


The McVittie metric provides a good approximation to 
the exterior metric of a massive spherical body with a 
physical radius much larger than its Schwarzschild radius. 
In the interior it is most convenient to work with 

R=[ 1 + fA (1 + jkr 2 )~ 1 / 2 a{t)r 

as the radial coordinate. If k = 0, then the surface 
{t, R) = const has area AttR 2 . When t = to the inte¬ 
rior geometry will be approximately Schwarzschild over 
scales where dl? ~ 0{R S ) and df ~ 0{Li). We define 
new dimensionless coordinates by 

T = {t — to)/Li and £ = R/R s . 

In ref. Ilfill we found the interior solution for the dilaton 
field to 0{5 2 ) and the exterior solution to order 0{S), 
and explicitly performed the matching. The inner ap¬ 
proximation to <f> is then given by: 

h-^ + s^P + o^), 

its 


where 



^ (5 (4 0) " + wf' - h'F (4)) + l^k{ST)F (0 O ) - h 2 F (fo) + (h'(ST) + k {ST )) (9) 

-\f (to) B M {ct> E ) RWfust + F (to) RlVu (0b)) + |r (^p- + D{Sp + O (|j) , 


4> c {t) = 1 i^E { T ) i s I^e cosmological value of dilaton field, 
and C{ST) and F (to) are ‘constants’ of integration to be 
specified by a boundary condition at R = R s depending 
on the details of the energy-distribution in R < R s . We 
do not expect these two latter terms to be any larger than 
the other terms in the above expression when £ = 1. Our 
primary concern is the behaviour of the time derivative 
of 4> in the interior and in the McVittie background. We 
have 


metric is 

ds 2 = df 2 - R ’ r ^^ dr 2 - R 2 (t, r){d9 2 + sin 2 Odcf 2 }, 
1 — k(r) 


where 


R 2 t = ~H r ) + 


2m 


t- 2 Z (r) 1 , ~ 

-- — + -A R 2 . 

R 3 


dt<t>i 

dttc 


1 - F ,0 {(j> 0 ) In 



Since, in most cases of interest, the dilaton to matter cou¬ 
pling is weak and we are far outside the Schwarzschild ra¬ 
dius of the Sun, both F((to) an d 2m/R are <C 1. Hence, 
dt4>i / dt4> c ~ 1 , and the local time evolution of the dilaton 
tracks the cosmological one. The strength of this result 
arises partly from our restrictive choice of background 
metric. 

We now consider Tolman-Bondi backgrounds, where 
more interesting behaviour is allowed. The background 


The matter content of these models is pressureless dust 
and there is a cosmological constant, A. To reproduce our 
required geometrical set-up, with the dust density posi¬ 
tive everywhere, and approach to spatial homogeneity at 
large r, we need 

Z{Rs)= 0, z, r >0, fim Z(r) = p d °l t H 2 r 3 . 

In the interior we take T = {t — to)/T/, and £ = R/R s 
as coordinates and 

Z(r)/m ~ S q ^i{r]) + o{S q ), 
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where 77 = (£ 3 / 2 + 3T/2) 2 ^ 3 and R s rj = r + 0(S q ,S 2 ^ 3 ). 
From the exact solutions (with r = R at t = to) we find 

k(rf) = 5 2 ' 3 k 0 (1 + 6 q n 1 ( 77 ) + o (S q )) + O (<5 5/3 ) , 

where fco = (2m/R s )(n / H 0 t) 2 / 3 . Solving for <f> in the in¬ 
terior to order 5 q we find 

<j> I ^(j>f ) (5T) + 8 q (j) { p +o(8 q ), 


where cjJf' 1 is given by 

4 0) = ^ (Vl - S 2 / 3 koT*^j - F (0o) in (l - ^ , 

and the transformation T —> T* is given in ref. [lbj 
along with the expression for 0^. We find that the time 
variation of <f>i is given by 


J 


^7,r, « <5^1 - S^ko^iVl - P/ 3 k 0 5T*) (l - Fj (4) In (l - 


( 10 ) 


-S q B^(t> c ) 


2m 

-Rs 


rv <w ^ 


1 Fit 0 )\ ^($ = 1 ,n 

^/2(v,t) 1 v s(0c)y c 3/2 


Here, the excluded terms are guaranteed to be smaller 
than the one of the two included terms; however, since 
they might be larger than the other term, both numeri¬ 
cally and in the limit S —> 0 , the above expression is not 
a formal asymptotic approximation. When the first term 
in eq. m dominates, condition © holds, and the local 
time evolution of the scalar field is, to leading order, the 
same as the cosmological one. Assuming our background 
choice is suitable for the application of this method, con¬ 
dition © fails to hold if, and only if, the second term in 
this expression dominates over the first. 

We have explicitly checked that the McVittie back¬ 
ground matching works, and by analysing the forms of 
the interior and exterior expansions in the Tolman-Bondi 
metric we can show that the matching procedure will 
break down if m < 1 —pa < 1—p(l—m) _1 or if m > 1 —pa 
when pa = 0, where pi (x) ~ x" as x ~ > °o for some 
n > 0, and where Z = Zfrw + b p z\ + 0(S q ), and 
Zi(H 0 r) ~ (H 0 r) m . By considering the behaviour of the 
next-to-leading order terms we can say that if both the 
interior and exterior zeroth-order approximations are si¬ 
multaneously valid in some intermediate scaling region 
we need 


lim (i? 2 Ke(i?, f)) = o(l) 

for all a £ (0,1); where the limit is taken with 

{L“L)j -a (t — t 0 ), LfL^~ a R} held fixed. This is also the 
condition for matching that we have found in the non- 
spherically symmetric extension of our result using the 
Szekeres spacetime 0- If such an intermediate region 
does indeed exist then the zeroth-order matching will be 
valid. The general form of the interior approximation to 
O ( S p ) will then be correct; the only unknown function 
in it is B(ST). If the matching works to order S p , as well 
as zeroth-order, then we have argued in |ljj| that B(ST ) 


will be quasi-static. If the matching procedure does not 
work to this order then its quasi-static character may 
be lost. However, we would not expect it to vary in time 
any faster than the other O (5 P ) terms. So long as we can 
match the zeroth-order approximations in some region, 
we can find the circumstances under which condition © 
holds by comparing the sizes of the two terms in W- 
By adding sub-leading order terms to this condition so 
as to express it in terms of physical quantities, we find 
fn| that to leading order it is equivalent to 


/ := 


(j) — <j) c B t< j,((f>c) Xy(ij) A(i? j (K;e) 
0c <p c (t) 



( 11 ) 


where dl = dr'I?, r , and the path of integration, "f(R), 
runs from R to 00 , and A(i?,te) = R,tne — HRne c (t ), 
A(i? t ) = R t — HR. R,t(r, t) is instantaneous the velocity 
of a matter particle at (r, t). Since we expect <f> will only 
feel its causal past, we should take 'y(R) to run from R 
to spatial infinity along a past radially-directed light-ray. 
Eq. (Hill has the form expected for condition © to hold 
inside a local, spherically-symmetric, inhomogeneous re¬ 
gion produced by a Schwarzschild mass embedded in an 
asymptotically FRW universe: the LHS of eq. GH> van¬ 
ishes as we approach the cosmological region and the lo¬ 
cal time variation of the dilaton held is driven by the time 
variation of the energy density along its past light-cone. 
The second, or drag, term in IIII comes from the change 
of time-coordinate T T* and is small compared to the 
first term whenever there is a local overdensity of matter. 

Therefore we are led to conjecture that eq. CCD is a 
general sufficient condition for eq. © to hold: applica¬ 
ble even in spherically-symmetric (or almost spherically- 
symmetric), dust plus cosmological constant, back- 
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grounds in which the matching procedure formally fails. 

If the cosmic evolution of the dilaton field is dominated 
by its coupling to matter, so \B^{$ c )ne c \ » \V phi {^ c )\, 
then <j> c (t) = O (B^((j) c )Ke c (t)/H{t)) and so condition © 
holds at the time t = to whenever 

/ dZff(io)!^^«l. ( 12 ) 

J'y(R) £ c\ t OJ 

where we assumed the drag term to be negligible and 
dropped it. Alternatively, if the cosmic evolution is 
potential-driven, so |V^(<(> C )| |B 0 (</> c )«e c |, then the 
LHS of the above expression will be suppressed by an 
additional factor of |.B 0 (</> c )Ke c /V^(^> c )| <C 1 . 

We can apply out results to solve our original prob¬ 
lem: to determine if condition © holds when the local 
region is the Earth or the solar system. Consider a star 
(and associated planetary system) inside a galaxy that is 
embedded in a large cluster. The cluster, of size R c iust , 
is assumed to have virialised. We consider the different 
contributions to the LHS of condition m in this astro¬ 
nomical hierarchy, with: I = I c i ust + I ga i + I star . 

Since we have only performed our calculations for 
spherically-symmetric backgrounds, we should, strictly 
speaking, also require spherical symmetry about our star. 
However, in a forthcoming paper we will show, by deter¬ 
mining the explicit matched asymptotic expansions for 
the case of the Szekeres background metric m , which 
possesses no Killing vectors [lf|, that criterion itTTl) still 
gives a very good approximation to the true condition 
in the absence of spherical symmetry. For illustration 
we have evaluated I for Brans-Dicke theory. For other 
theories, I will still take very similar values. We assume 
that the density of non-virialised matter, just outside the 
virialised cluster, is no greater than the average density 
of the cluster and if we move away from the edge of the 
virialised region, the over-density drops off quickly, i.e. 
as R~ s , s > 1/2. In this case, the Tolman-Bondi result, 
eq. US. applies and the magnitude of I is bounded by: 

hlust < to 1 (S - 1/2)"V2 MdustRdust^ 1 

£c 

^ 1.2 X 10 I R c lust \ / Vclust \ / | ,3 

1/2) \ Mpc / Vkms -1 / ™ 

» 1.5 x 10 _ 3 [ 3 /( 2 s - l)]n- x (l + z vir f < 1 , 

where we have used 3 GM/R c i ust = v 2 lust ; R c iust is the 
post-virialisation scale of the cluster, and v c iust is its av¬ 
erage virial velocity; to = li.lGyrs is the age of the 
universe. In the final approximation we used represen¬ 
tative values Rciust = 100 Mpc and v c i us t = 200 km s ^ 1 
appropriate for a cluster like Coma. Taking a cosmo¬ 
logical density parameter equal to fl m = 0.27, in accor¬ 
dance with WMAP we expect that for a typical clus¬ 
ter which virialised at a redshift z V i r <C 1 , we would 
have I c iust ~ 5.7 x 10~ 3 . The term in [] is unity when 
s = 2, i.e. 2 GM/r —> const ; such a matter distribution 


is characteristic of dark matter halos. Different choices 
of s > 1/2 can be seen to only change this estimate by 
an 0(1) factor. As s —> 1/2 the matched asymptotic 
expansion method, and hence this particular evaluation, 
breaks down. We believe that the generalised formula 
for I will, however, still give accurate results. For the 
star and galaxy contribution, we assume each is accret¬ 
ing matter from the interstellar medium (ISM) and in- 
tergalactic medium (IGM) respectively. We find 

Igal « — 2tg 1 2GMgaiR ( Q al ' > 

Mq \ v clust ) 

« —4.4 X 10~ 6 (1 + Zyir ) 3 1, 


<-)« — 1 /r)/^i ^-ISM 

— 2 1 0 \J2LxmrG - 

£c 

-1.2 xlO" 13 !!- 1 -^ 

m Mq 

-1.8 x 10” 14 n< 1, 


( kms 1 
VlSM 


h 2 n 


where we have taken cigm ~ Vclust and vigm ~ 
Vclust, Mgai = 10 12 Mq, Vclust = 200kms -1 , R G = 
2 GM ga i/v 2 IGM , v g = 2 Gm/v 2 ISM , and cism = n 
protons cm -3 , where we expect n ~ 1 — 10 4 , using 
vism = 5km s^ 1 and m = Mq for the numerical es¬ 
timates. We take h = (Ho/100kms~ 1 Mpc~ 1 ) = 0.71 
for the Hubble constant. It is clear that, in general, 
1 1 star I < 14a/1 < |4/u«i| as expected. The infall of dust 
into the cluster will tend to be the dominant contribu¬ 
tion to the LHS of eq. f*H?l) . Our estimate of I c i U st should 
be viewed as an upper bound on its value, which is very 
small compared with unity and hence we expect condi¬ 
tion 0 to hold near the Sun; with deviations from that 
behaviour bounded by the 0.6% level (if the dilaton cou¬ 
ples only to baryonic matter this could rise to as high 
as 3%). Assuming that the conditions in our solar sys¬ 
tem are not too different from those considered above, we 
conclude that irrespective of the value of the dilaton-to- 
matter coupling, and what dominates the cosmic dilaton 
evolution, that the condition © will hold in the solar 
system in general, and on Earth in particular, to high 
precision. Cosmological variations of ’constants’, like a, 
will therefore be detected directly in terrestrial experi¬ 
ments to the accuracy specified by the value of I. 

We can also conclude that, even if condition © is vi¬ 
olated during the collapse of an overdense region of mat¬ 
ter, then once the region stops collapsing as a result of its 
virialisation, the time evolution of the dilaton field con¬ 
tinue evolving towards homogeneity. Contrary to what 
has been assumed before in the literature, local virial¬ 
isation does not stabilise the value of the dilaton, and 
protect it from any global cosmological variation. In ad¬ 
dition, by considering what happens when 2 m/R s —> 0, 
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we show in ref. |l6j ] that we expect the time variation 
of </> (and hence G in a scalar-tensor gravity theory like 
Brans-Dicke ,19] ) on the horizon of a black-hole (as mea¬ 
sured by an infalling observer) to track its cosmologi¬ 
cal value, confirming the prediction made by Jacobson, 
|20j . the inhomogeneous models of ref. |2J, and the 
numerical calculations of Harada et. al., [9|, who also 
studied the Tolman-Bondi background. We have found 
that whenever weak condition, e;/e c <C l/(2mifo)(i§> 1) 
holds, that Jacobson’s particular solution to Ocj> = 0 in 
the Schwarzscliild metric is a true asymptotic approxi¬ 
mation to the local behaviour of the dilaton close to the 
Schwarzschild horizon. 

In conclusion: we have used the method of matched 
asymptotic expansions to find a sufficient condition for 
the local time-variation of a scalar field, or varying phys¬ 
ical ‘constant’ driven by such a field, to track the cos¬ 
mological evolution, and we have proposed a generali¬ 
sation of this condition that is applicable to scenarios 
more general than those explicitly considered here. By 
numerically evaluating the condition for an astronomical 
scenario similar to the one appropriate for our solar sys¬ 
tem, we have concluded that almost irrespective of the 
form of the dilaton-to-matter, and the form of the dila¬ 
ton self-interaction, its time variation in the solar system 
and in terrestrial laboratories will track the cosmolog¬ 
ical one. We have therefore provided a general proof 
of what was previously merely assumed: that terrestrial 
and solar system based observations can legitimately be 
used to constrain the cosmological time variation of many 
supposed ‘constants’ of Nature. 
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